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In this paper we study the asymptotic behaviour of power and multipower variations of pro- 
cesses Y: 



Y t = [ g(t-s)a s W(ds) + Z t , 

J — oo 

where g : (0, oo) — > R is deterministic, a > is a random process, W is the stochastic Wiener 
measure and Z is a stochastic process in the nature of a drift term. Processes of this type serve, 
in particular, to model data of velocity increments of a fluid in a turbulence regime with spot 
intermittency a. The purpose of this paper is to determine the probabilistic limit behaviour of 
the (multi)power variations of Y as a basis for studying properties of the intermittency process a. 
Notably the processes Y are in general not of the semimartingale kind and the established theory 
of multipower variation for semimartingales does not suffice for deriving the limit properties. 
As a key tool for the results, a general central limit theorem for triangular Gaussian schemes is 
formulated and proved. Examples and an application to the realised variance ratio are given. 

Keywords: central limit theorem; Gaussian processes; intermittency; non-semimartingales; 
turbulence; volatility; Wiener chaos 

1. Introduction 

The motivation for the development of the results reported in this paper has been the need 
to construct tools for studying the probabilistic limit behaviour of (realised) quadratic 
variation and other multipower variations in relation to the class of Brownian semistation- 
ary (BSS) processes. This class, which was introduced in [12], consists of the processes 
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Y = {Y t }t£R that are defined by 



Y t = n + 



I 



— OO 



g(t - s)<r s W(ds) + 



— OO 



q(t — s)a s ds 



(i.i) 



where /i is a constant; W is a Brownian measure on R; g and q are non-negative deter- 
ministic functions on K, with g(t) = q(t) = for t < 0; and a and a are cadlag processes. 
When a and a are stationary, so is Y. Hence the name Brownian scmistationary processes. 
It is interesting to note that the fractional Ornstein-Uhlenbeck process is, in fact, also 
representable in the form (1.1). The same is true of a wide class of stable pseudo- moving 
average processes; see Corollary 4.3 in [2]. 

The BSS processes form the natural analogue, for stationarity related processes, to 
the class BSM. of Brownian semimartingales 



In the context of stochastic modelling in finance and in turbulence, the process a em- 
bodies the volatility or intermittency of the dynamics, whether the framework is that of 
BSAi or BSS. For detailed discussion of BSS and the more general concept of tempo- 
spatial ambit processes see [7-12]. Such processes are, in particular, able to reproduce 
key stylized features of turbulent data. 

A main difference between BSAi and BSS is that, in general, models of the BSS form 
are not semimartingales (for a discussion of this, see Section 3 of [12]). In consequence, 
various important techniques developed for semimartingales, such as the calculation of 
quadratic variation by Ito calculus and those of multipower variation, do not apply or 
suffice in BSS settings. The present paper addresses some of the issues that this raises. 

The theory of multipower variation was primarily developed as a basis for inference 
on a under BSA4 models and, more generally, Ito semimartingales, with particular focus 
on inference about the integrated squared volatility a 2+ given by 



This quantity is likewise a focal point for the results discussed in the following. 

Section 2 introduces common notation for multipower variation and recalls some basic 
properties of such quantities. A law of large numbers and a central limit theorem for 
multipower variation of triangular arrays of Gaussian random variables are derived in 
Section 3, and these limit results are drawn upon in Section 4 to establish probability 
and central limit theorems for multipower variation for BSS processes, with most of 
the proofs postponed to Section 8. Section 5 presents several examples and Section 6 
discusses an application concerning the limit behaviour of the realised variation ratio, 
that is, the ratio of realised bipower variation to realised quadratic variation. Section 7 
concludes and indicates some possible directions for further related work. 




(1.2) 




(1.3) 
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2. Multipower variation 

The concept of (realised) multipower variation was originally introduced in [13] in the 
context of semimartingales, and the mathematical theory has been studied in a number 
of papers [6, 17, 24, 27] while various applications are the main subjects in [14-16, 25, 33]. 
Multipower variation turns out to be useful for analysing properties of parts of a process 
that are not directly observable. In this section we present the definition of realised 
multipower variation and recall its asymptotic properties for some classes of processes. 

Let us consider a continuous-time process A, defined on some filtered probability space 
(O, J 7 , (j r t)t>0j P), that is observed at equidistant time points U = i/n, i = 0, . . . , [nt]. A 
realised multipower variation of the process A is an object of the type 



l\K +J -iX\ P3 , ^X = X i/n -X {i _ x)/n , Pl ,...,p k >0, (2.1) 



for some fixed number k > 1. We now present an overview of the asymptotic theory for 
quantities of the form (2.1) for various types of processes X. 
We start with the BSM case 



where W is a Brownian motion, a is a locally bounded and predictable drift process and 
a is an adapted and cadlag volatility process. As was established in [6], the convergence 
in probability 



[nt]-fc+l k 



3=1 




(2.2) 




(2.3) 






(2.4) 
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where B is another Brownian motion, defined on an extension of the probability space 
(CI, J-, (J~t)t>o, P) and independent of J 7 , and the constant A is given by 

k k k—1 rn k k—m 

A=Y[(jQ Pl -(2k-i)Y[i4 l + 2 Yl Tl^pi n n ^pi+pi+m- 

1=1 1=1 m=l/=l l=k-m+l 1 = 1 

Recall that the stable convergence of processes is defined as follows. A sequence of 
processes Z n converges stably in law towards the process Z (written Z n — > Z), that is, 
defined on the extension of the original probability space (Q, J 7 , (J r t)t>o, P), if and only if 
for any bounded and continuous real-valued functional / and any ^-measurable random 
variable V it holds that 

lim E[f(Z n )V]=E[f(Z)V], 

n— f oo 

we use the notation Z n —> Z. 

A crucial property of the realised multipower variation is its robustness to jumps when 
ma,Xi(pi) < 2 [17, 26]. Assume for a moment that A is a general Ito semimartingalc 
with continuous part X c satisfying (2.2). Then, by (2.3) and the robustness property, we 
obtain the convergence 

^ 2E M-l |A » x||A „ +iX| ^ [xc] 

EK|A?x| 2 W 

where [X] denotes the quadratic variation of the semimartingalc X and the limit is less 
than or equal to 1. The latter result, together with the stable convergence in (2.4), can 
be used to construct a formal test for jumps (see [13]). On the other hand, we know that 
if the limit of the left-hand side is greater than 1 (which is the case for some typical 
turbulence data), the process X cannot be an Ito semimartingalc. 

In another direction, a study [5] was made of the asymptotic behaviour of bipower 
variation for processes of the type 

A t = A + / a s dG s , t>0, (2.5) 
Jo 

where G is a continuous Gaussian process with centered and stationary increments (the 
latter integral is defined as a Riemann-Stieltjes integral). The process defined in (2.5) is, 
in general, also not a semimartingale, and the theory in [6] does not apply. In particular, 
a different normalisation is required. Define the (normalised) multipower variation by 

[nt]-fc+l k 

V(X,pi, . . . ,Pk)t = — E I|A? +3 -_iX|w, Pi,..-,Pk>0, 



where t„ > is given by 



rl=R[{l/n) 



(2.6) 
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with 

R(t)=E[(G s+t -G s ) 2 }. (2.7) 
Under some assumptions on R and the volatility process a it was shown that 

V(X, Pl ,...,p k )?^p Pu ... ;Pk /VsPds 

Jo 

for a certain constant p Pll ... lPk that depends on the behaviour of R near 0. Furthermore, 
an associated (stable) central limit theorem, of a form similar to (2.4), was derived. Note, 
however, that in general there are essential differences between the characters of BSS 
processes and processes of type (2.5). In the latter case, the process a has only a local 
influence in the value of X whereas, in the BSS case, the process is also affected by the 
past of a. 



3. Multipower variation of Gaussian triangular arrays 

In this section we derive some asymptotic results for functionals of arrays of stationary 
Gaussian sequences. We consider a triangular array (-Xi,n)n>i,i<i<[nt] (t > 0) of row-wise 
stationary Gaussian variables with mean and variance 1. Let 

r n (j) = cor(Xi i „,Xi +J> ), j > 0, (3.1) 

be the correlation function of (^Q,n)i<i<[nt] • Assume that the array (^i.n) n >i,i<j<[nt] 
is "non-degenerate", that is, the covariance matrix of (Xi iTl , . . . ,Xi +k ,n) is invcrtible for 
any k > 1 and n > 1 (otherwise the results below do not hold) . 

Now, define the multipower variation associated with the sequence (-Xi,n)n>i,i<i<[nt] : 

[nt]-k+l k 

V(p u ...,p k ) n t =- ]T IJ|X i+i _ lin |«, Pi,---,Pk>0- (3-2) 

H »=1 3=1 

Our first result is the weak law of large numbers. 
Theorem 1. Assume that there exists a sequence r(j) with 

n— 1 

r£(j)<r(j), -^ r (j)^0 (3.3) 

3=1 

as n — > oo . T/ien it holds that 

V{px,... t p k yi-p^] ..., Pk t^0, (3.4) 

where 

P P t.., Pk =E[\Xi, n \ Pl ---\Xk,n\ Pk }. (3.5) 
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Proof. See Section 8. □ 

Before we present the associated central limit theorem, we need to introduce another 
Gaussian process. Suppose that r n (j) — > p(j), j = 1, . . . , k — 1, for some numbers p(j). 
Let (Qi)i>i be a non-degenerate, stationary, centered (discrete-time) Gaussian process 
with variance 1 and correlation function 

p(J) = cat(Qi,Q 1+j ), j>l. (3.6) 

Define 



V Q (p 1 ,...,p k )? = ± J2 fl\Qi+i-i\ Pi (3-7) 

»=1 j=l 

and let p Pl ,... tPk = E{\Q 1 \Pi---\Q k \P«}. Then ppi}...,p k pp 1 ,...,p k and in this case we obtain 
the uniform convergence on compacts in probability (or ucp convergence): 

V(pi,...,Pk) t — > p Pu .... Pk t- 

Now we have the following central limit theorem for the family (V{p\, . . . ,Pk)f)l<j<d- 

Theorem 2. Assume that 

r»(j)->p(j), J>0, (3.8) 
and that, for any j,n> 1, there exists a sequence r(j) with 

oo 

r 2 n (j)<r(j), £r(j)<oo. (3.9) 

Then we have 

V^(V(A,---,rt)?-p ( f v ,t)i< 3 < d ^P 1/2 Bu (3.10) 
where B is a d- dimensional Brownian, j3 is a d x d-dimensional matrix given by 

Pa= lim ncov(VQ(pi,...,piK,VQ(p> v ...,p j k ) 7 !), l<i,j<d, (3.11) 

n— f oo 

and the weak convergence holds in the space 2?([0,T] d ) equipped with the uniform topology. 
Proof. See Section 8. □ 
Remark 1 . Similar asymptotic results can be obtained for general quantities of the form 

^ [nt]-fc+l 

— y H(Xi n , . . . ,Xi + k-l.n) (3-12) 
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for some function H : R fc — > R. Let m denote the Hermite index of H (notice that the 
Hcrmitc index of the power function used in (3.2) is 2). Replace condition (3.3) by 

n— 1 

IC(i)|<r(i), -£r(j)^0 
i=i 

and (3.9) by 

IC(j)l<r(j), ]>>(j)<oo. 

j=i 

Then Theorems 1 and 2 hold true for the functional (3.12) provided that EH 2 (N k (0, £)) < 
oo for any invertible £ <E R fcxfc . We omit the details. 

Remark 2. Ho and Sun [23] have shown a non-functional version of Theorem 2 for 
statistics of the type (3.12) when the correlation function r n docs not depend on n. To 
the best of our knowledge, Theorem 2 is the first central limit theorem for (general) 
multipower variation of a row-wise stationary Gaussian process. 

4. Multipower variation for BSS processes 

Armed with the general theorems proved in Section 3, we are now set to establish laws 
of large numbers and central limit results for multipower variations in the framework of 
the Brownian scmistationary processes. The regularity conditions invoked are given in a 
first subsection, while the next states the theorems, the main parts of the proofs being 
postponed to Section 8; the third subsection discusses the nature of the rather technical 
regularity conditions and describes a set of simpler assumptions that are more amenable 
to checking. 

4.1. Conditions 

We consider a filtered probability space (f2, J 7 , (J-t)t>o, P), assuming the existence thereon 
of a BSS process, without drift term for the time being, that is, 

Y t = [ g(t-s)a s W(ds), (4.1) 

J — oo 

where W is an (J r t )-Brownian measure on R, a is an (J-" t )-adapted and cadlag volatility 
process and g : M. — > R is a deterministic continuous memory function with g{t) = for 
t < and g £ £ 2 ((0, oo)). We also require /_ g 2 (t — s)a 2 s As < oo a.s. to ensure that Y t < 
oo a.s. for all t > 0. By an (J r i)-Brownian measure we understand a Gaussian stochastic 
measure such that, for any Borelian set A with E^W^A) 2 ] < oo, 



W(A)~N(0,m{A)), 
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where m is the Lebesgue measure, and if A C [t, +00), then W(A) is independent of Ft- 
Note that {B t := f W(ds),t > a} is a standard Brownian motion starting in a. 

The process Y is assumed to be observed at time points ti=i/n, i = 1, . . ., [nt]. Now, 
let G be the stationary Gaussian process defined as 

G t = f g(t-s)W(ds). (4.2) 
J — 00 

This is an important auxiliary object in the study of BSS processes. Note that G belongs 
to the type of processes occurring in (2.5), and that the autocorrelation function of G is 

r(t) = Jo y r l 2 /r, ■ (4-3) 

Jo r(«)a« 

We are interested in the asymptotic behaviour of the functionals 

[nt]-fe+l fe 

" T " i=l i=l 

where A"F = F i/n - and = R(l/n) with i?(t) = E[\G s+t - G s \ 2 }, t > 0. In the 

following we assume that the function g is continuously diffcrcntiablc on (0,oo), \g'\ is 
non-increasing on (6,00) for some b > and g' € L 2 ((e,oo)) for any e > 0. Moreover, we 
assume that for any t > 0, 

/oo 
(ff , (a))V t a _.da<oo (4.4) 

almost surely. 



Remark 3. Assumption (4.4) ensures that the process Y has the same "smoothness" as 
the process G (see Lemma 1 in Section 8). It is rather easy to check in practice, because 
it is implied by the condition EF t < 00 for t > 0. Furthermore, if g has bounded support, 
assumption (4.4) is trivially fulfilled since a is cadlag. 



Remark 4- Let us note again that the process Y is, in general, not a semimartingale. 
In particular, this is the case when g' ^ L 2 ((0,oo)). For a closer discussion, see [12]. On 
the other hand, the process Y is not of the form (2.5). Thus, we require new methods to 
prove the asymptotic results for V(Y,pi, . . . ,Pk)t ■ Processes of the form (4.1) are used 
for modelling velocity of turbulent flows; see [8-10]. The function g, which is used in 
such models, behaves often as x s near the origin. Hence, when S £ (—1/2, l/2)\{0}, Y is 
neither a diffcrcntiablc process nor a semimartingale (because g' £ L 2 ((0,oo))). This is 
the primary case of our interest. 
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We define the correlation function of the increments of G: 

rnU) = covf^, A " +jG ) = m + l)M + m-DM-2RU/n) t > Q 

Next, we introduce a class of measures that is crucial for our purposes. We define (recall 
that g(x) := for x < 0) 

?r"(A) = ^ 9(X - 1/n) - g(x))2 dX , A e B(R). (4.5) 
S™{9{x- 1 M -g(x)) 2 dx' 

Note that 7r" is a probability measure on R + . 

For the weak law of large numbers we require the following assumptions: 

(LLN) There exists a sequence r(j) with 

^ n— 1 

r*(j')<r(j), -^r(j)^0. 
n * — ' 

i=i 

Moreover, it holds that 

lim 7r"((£,oo)) = (4.6) 

n— »oo 

for any e > 0. 



Remark 5. 

(i) The first condition of (LLN) is adapted from Theorem 1. It guarantees the ucp 
convergence of V(G,pi, . . . ,Ph)t- The second condition of (LLN) says that the 
whole mass of the measure ir n concentrates at 0. In particular, it is equivalent to 
the weak convergence 

where 5q is the Dirac measure at 0. 

(ii) Condition (4.6) is absolutely crucial for the limit theorems given in the next sub- 
section. When this condition is violated things become more complicated. In par- 
ticular, it may lead to a different stochastic limit of V(Y,pi, . . . ,pk)'t ( see the first 
example in Section 5). Intuitively, this can be explained by the observation that 
the increments A™Y contain substantial information about the volatility (far) out- 
side of the interval [^-, ^] when condition (4.6) does not hold. Thus, in general, 
we can not expect the limit described in Theorem 3 below. 



Now, we introduce the assumptions for the central limit theorem: 
(CLT) Assumption (LLN) holds, and 



r n (j)->p(j), 



3>0, 
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where p(j) is the correlation function of (Qi)i>i, as introduced in (3.6). Furthermore, 
there exists a sequence r(j) such that, for any j, n > 1, 

oo 

rl(j)<r(j), E r (i)<°°, 
i=l 

and for some 7 € (0, 1] we have 

E{\a t -a s \ A }<C\t-s\ A ~< (4.7) 

for any A > 0. Finally, set p = mini<i<fc i<j<d(p^). Assume that j(p A 1) > \ and that 
there exists a constant A < — such that for any e n = 0(n,~ K ), k G (0, 1), we have 

7r"(( £ „,oo)) = 0(n A ( 1 - K )). (4.8) 



4.2. Limit theorems 

In this section wc present the limit laws of multipower variations of BSS processes, 
in part widening the scope slightly to allow more general drift terms. Recall that the 
(realised) multipower variation of a process Y of the form (4.1) is defined as 

[nt]-k+l k 

V(YiPi,...,Pk)? = — p+ E Il A ^-i^r. Pi,~-,Pk>0, (4.9) 

t=l 3=1 

where t% — R(^) and p + = X^=iPj- Our first result is the following probability limit 
theorem. 



Theorem 3. Consider a process Z = Z\ + Z2, where Z2 — Y is given by (4-1)- Assume 
that the condition (LLN) holds and that 



[nt]-k+l k 

^ e n^v^r^o, 

ft' In 1 1 

i=l j = l 

where i\, . . . , t& € {1, 2}, /or any t > and any (ti, . . . , t&) ^ (2, . . . , 2). Define 



(4.10) 







Pi 

















Then we have 



V(Z. 



Jo 



ucp 



Proof. Sec Section 8. 



(4.11) 

□ 



Asymptotic results 



1169 



Remark 6. The multipower variation is robust to drift processes Z\ that are smoother 
than the process Y. Assume, for instance, that the process Z\ satisfies 

E[\Z 1 (t)-Z 1 (s)\?}=o(Ri>/ 2 (\t- S \)) 

for every p > 0. In this case, condition (4.10) is obviously satisfied. 

Next, we demonstrate a joint central limit theorem for a family [V{Z 1 p\, . . . ,P J k )t)i<j<d 
of multipower variations. 

Theorem 4. Consider a process Z = Z\ + Z2, where Zi = Y is given by (4-1)- Assume 
that the condition (CLT) holds and that 

[nt]-k+l k 
y/nT n + i=l 1 = 1 

where o±,...,if. £ {1, 2} ; for any t > and any (ti, . . ., tfc) 7^ (2, . . . , 2). Then we obtain 
the stable convergence 

^(v(Z,p{,...,pi)?-p% ) fla^ds) -%f A y^ Bs , (4.12) 

V Pi>-'PkJ / l<j<d Jo 

where B is a d-dimensional Brownian motion that is defined on an extension of the filtered 
probability space (Q,,F,(Ft)t>o,P) an d is independent of T , A is a dx d-dimensional 
process given by 

A«=A i |<7 8 |fW +) l<i,j<d, (4.13) 
and the dx d matrix ft is defined in (3.11). 

Proof. See Section 8. □ 
4.3. Discussion of assumptions 

We start our discussion again by considering the auxiliary, centered, stationary Gaussian, 
process 

G t = / g(t-s)W(ds). 
J —00 

First of all, we want to demonstrate how Theorems 1 and 2 apply for the multipower 
variation of the process G. In other words, we will give a hint how to check the conditions 
of these theorems. 
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Recall definition (2.7) of the variance function R of the increments of G and note that 
R(t)=E[\G s+t -G s \ 2 } = / g 2 {x)dx+ / (g(t + x) - g(x)f dx, *>0. 



Clearly, the asymptotic behaviour of the multipowcr variation of the process G is fully 
determined by the behaviour of the function R near 0. As we deal with a continuous 
process G, it is natural to assume that R(t) behaves essentially as t a (for some a > 0) 
near (later on we will formalize this assumption). Since the case where the paths of 
G are differentiablc (a.s.) is not very interesting for us (because the consistency can 
be deduced by the mean value theorem), we concentrate on the region < a < 2 (the 
corresponding g(t) behaving as ^ a_1 ^' 2 ). 

Let us introduce a new set of assumptions that correspond to the previous discussion. 
These assumptions were proposed by Guyon and Leon in [22] (those authors considered 
the case of centered, stationary Gaussian processes X; this relates to the BSS setting 
with a constant) and the same assumptions were used in [3] and [5]. 

(Al) R(t) = t a Lo(t) for some a £ (0, 2) and some positive slowly varying (at 0) function 
I/O) which is continuous on (0,oo). 

(A2) R"(t) = t a ~ 2 Li(£) for some slowly varying function L2, which is continuous on 
(0,oo). 

(A3) There exists b £ (0, 1) with 



K = lim sup sup 

x~¥0 y£[x,x b ] 



< 00. 



L (x) 

Recall that a function L : (0, 00) — > M. is called slowly varying at when the identity 

hm — — = 1 

x\a L{x) 

holds for any fixed t > 0. 

Now, note that under assumption (Al) we have, for any j > 1, 



r n {j) = cov 



= R((J + 1)/") + R((J - 1)/") - 2RUM (j) u 14) 

2R(l/n) P[J) [ ' ' 

= \((j+ir-2r+(j-ir), 

because Lq is slowly varying at 0. It is obvious that p(j) is the correlation function of the 
discrete-time stationary Gaussian process Qi = B^ 2 — Bf^ 2 , where B a l 2 is a fractional 
Brownian motion with parameter a/2. 
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Remark 7. It is easy to see that the convergence of 

R((j + I) /n) + R((j - l)/n) - 2R(j/n) 
2R(l/n) 

to some real number for all j > 1 implies that there exits an s(j) such that 

Since the result in Theorem 4 is independent of the scale of time we use, we must have 



R(A) aid 



s(j) 



for any A and then s(jk) = s(j)s(k); consequently s(j) = j a , for a certain a € R. More- 
over, since (j + l) a — 2j a + (j — l) a is a covariancc function, we have < a < 2. So in the 
present setting, (Qi)i>i, as defined in Section 3, is always a standard fractional Gaussian 
noise. 



As shown in [3, 22] assumptions (A1)-(A3) imply that condition (3.3) holds for any 
a £ (0,2) and condition (3.9) holds for any a € (0,3/2). Hence, Theorem 1 holds for all 
a e (0, 2) while Theorem 2 only holds for a € (0, 3/2). 

Now, let us see what the conditions (A1)-(A3) mean for the memory function g. For 
simplicity, let us consider functions of the form 

g(x) = x s l m {x), x>0. (4.15) 

For such functions we readily obtain assumptions (Al) and (A2) with 

a = 26 + l, S€ (-i,0)U(0,|) 

(the technical assumption (A3) has to be checked separately; for an example, see Sec- 
tion 5). Note that for 8 — 0, for which assumption (A2) does not hold, the process G is 
a semimartingale and the multipower variations can be treated as in [6]. 

Next, we discuss the assumptions of Section 4.1 for the function g defined in (4.15). 
Recall that condition (4.4) is automatically satisfied for functions g with compact support 
(as in (4.15)). A straightforward calculation shows that 

^((e,oo)) = 0((n £ ) 25 - 1 ) 

for any e > -k Thus, condition (4.6) of (LLN) is satisfied (because 26 — 1 < 0) and 
Theorem 3 is"valid for all 5 £ (-1/2,0) U (0, 1/2). 

Finally, we explain how to verify condition (4.8) of (CLT). Recall that p = 
mmi<i<k,i<j<d(Pi)- Let e„ — n~ K , k£ (0,1). We readily deduce that 



7r n ((e,oo)) = 0(n A ( 1 -' i )), X = 2S-l. 
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Thus, condition (4.8) is satisfied if 



1 Ap 

We immediately deduce that Theorem 4 holds if 

p>1: i>\, *ef-io), 

Remark 8. Clearly, we can deal with a larger class of functions g than g(x) = x 5 l(o.i] ( x ) ■ 
Assume that condition (4.4) holds. In the following we consider functions L g , L g >, which 
are continuous on (0,oo) and slowly varying at 0. We assume the following conditions: 
Assumption: g £ L 2 ((0,oo)) and for some 5g (-1/2, 0) U (0,1/2) it holds that: 

(i) g(x)=x s Lg(x). 

(ii) g'{x) = x s 1 L g i(x) and, for any e > 0, g' G L 2 ((e,oo)). Moreover, \g'\ is non- 
increasing on (b, oo) for some b > 0. 

We further assume that the function 

f-t poo 

R{t)= / g 2 (x)dx+ / (g(t + x) - g(x)f dx 
Jo Jo 

satisfies conditions (A1)-(A3) with a = 26 + 1. 

Under these assumptions we conclude (as for the simple example g{x) = a^l^i] (x)) 
that Theorem 3 holds for any S £ (-1/2, 0) U (0, 1/2), and Theorem 4 holds when further 

P>1: 7>i ^(-io), 
In both cases we have Qi = R 5+1 ^ 2 — B^V" 1 , i > 1. 



5. Examples 

This section discusses two examples of choice of the damping function g and the associ- 
ated probabilistic limit behaviour. 

As above let r denote the autocorrelation function of G = /_ g(- — s)W(ds). Note 
that assumptions (A1)-(A3) could equivalently have been formulated in terms of 1 — r 
rather than R (since R{t) = 2\\g\\ 2 (l - r(t))). 
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Suppose first that 



<?(*) 



-At 



l(o,i)(i) 



with A > 0. This example (for a detailed discussion, see [12]) is a non-semimartingale 
case, and it can be shown that 7r™ — > 7r, with n given by 



1 



7T 



-So 



1 



,2,\ 



Si, 



1 + c~ 2A " ' l + i 
where Si is the Dirac measure at i. Moreover, 

P r- , -2A\-1 2+ | ^ | ^A-i-l^l- 



V(Y,2)? 
where for any a < b 



(1 + e— )-Vr + (l + e-)-Vf+ M _ 1]5 



r 2+ 

>,6] 



cr ds. 



Thus, in particular, we do not have V(Y, 2)™ — > cr 2+ . Note that in this example assump- 
tion (A2) is not satisfied. 
Our main example is 

g(t)=e- 1 e- xt l {0tOo) (t) (5.1) 

for A > and with v > \. (So, for t near 0, g(t) behaves as t s with 6 = v — 1). The 
autocorrelation function is given by 



\2v-l 



r( t ) M. e"^ 

1 j ~r(2^-i) e 



(t + u) 



i/—l,,v—l n — 2 Xu 



du. 



(5.2) 



It can be proved using properties of Bcssel functions (see Sections 5.1-5.3 in [4]) that 
assumptions (A1)-(A3) are met provided that a = 2v — 1 G (0,2), that is, v 6 



'1 3> 



(n) 

and that p , , may be substituted by p j j in the central limit theorem provided 

r p{,...,pf, J u r Pi,—,p k 

v E ( i , | ) . For t->0 we also have the following asymptotic equivalence 



l-r(t). 



2 r>Ti72) 4 

|(Ai) 2 |logt| 
14(^)^ 2 + ^ 1 ) 



0(t 2 



* 1 3 
tor - < v < - , 

2 2 

tor t/ = — , 

2' 

- 3 

for - < v. 



Remark 9. So for ^ < j/ < 2 the autocorrelation function is twice diffcrcntiablc at and 
consequently Y" has continuously differentiable sample paths, while for ^ < ^ < | the 
sample paths are Lipschitz of order A for every < A < v — \ (cf. [21], Section 9.2). 
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6. An application 

Let us consider the realised variation ratio (RVR) defined for a stochastic process X as 

RVR?:=WPX±^. (6.1) 

The RVR is of interest as a diagnostic tool concerning the nature of empirical processes. 

In particular, it can be used to test the hypothesis that such a process is a Brownian 
scmimartingale (with a non-trivial local martingale component) against the possibility 
that it is of this type plus a jump process, see [17, 26] (some related work is discussed in 
[34]). If a jump component is present, then the limit of RVR" is smaller than 1. 

However, in the course of the turbulence project, mentioned earlier, when calculating 
the RVR for an extensive high-quality data set from atmospheric turbulence it turned 
out that the values of RVR were consistently higher than 1 . The wish to understand this 
phenomenon has been a strong motivation for the theoretical developments described in 
this paper. As a consequence of Theorem 3, we obtain the following probability limit 
result for the realised variation ratio of BSS processes: 

RVR?-i)(r n (l))^0, (6.2) 

where 



ip(p) = yl — p 2 + parcsinp, (6.3) 

which equals ^E{|C/V^|} of two standard normal variables U and V with correlation p. 
Moreover, we have that 

y/n(RVR t - ^(r n (l))) = Vn - r — ^ 

V Jo °s ds 

_ my, 2,0)?- rVd« 
- V^RVR? t — — — 

V Jo ds 

so, if the parameter a € (0, 1), by applying Theorems 3 and 4, we obtain 

MRVR?-TP(r n (l))) A (* ^(p(i))W4?Mr' ^ 6 - 4 ) 

\ ' Jo a s 

where ip is as above and the matrix /3 is given in Theorem 2. Specifically, we find (3 = 
Wij)i<i,j<2, where 

Pu= lim nvar(l/ Q (l,l)™), 

n— yoo 

I3 22 = lim nvar(y Q (2,0)5 l ), 
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0ia= lim ncov(y Q (2,0K,Vq(l,lK) 

n— >oc 

with Q as defined in Theorem 2. Thus, we obtain 

oo oo 

P 22 = var(Q?) + 2^ cov(Q?, Q\ +k ) = 2 + 4^ p 2 (k). 

k=l k=l 

Similarly, we have that 

oo 

/ 3 12 = 2cov(Q2 j |g 1 ||Q 2 |) + 2 ^cov(Q2 j |g 1+fc ||g 2+fc |). 

fc=i 

Then, if we write £[1X^X2X3!] := h(pi 2 , P13, P23), where Xi,X 2 and X3 are standard 
normal with cov(Xj,Xj) = pij, we have 

/8i a = ^(l,p(l),p(l))-^(p(l))) +2^^(p(*),p(fc + l),p(l))-|^(p(l))). 
To compute the latter, we may use the following formula: 

MPi2,Pi3,P23) = -(a/ 1 - Ph( 1 + P12 + P13) + (P23 + 2pi2Pi 3 )arcsin(p 2 3)); 

sec [28]. For the remaining term we deduce 

00 

/3n =var(|Q 1 ||Q 2 [) + 2^co V (|Q 1 ||Q 2 |,|Q 1+fc ||<g 2+fc |). 
fe=i 

However, while there is no explicit formula available for the latter expression, it can be 
computed numerically. 

7. Conclusion and outlook 

In this paper we have derived convergence in probability and normal asymptotic limit 
results for multipower variations of processes Y that, up to a drift-like term, has the form 

Y t = [ g(t-s)a s W(ds). 



A key type of example has g(t) behaving as t s for t \, and 8 e (—5, 2 1 )\{0}- I n those 
instances, Y is not a semimartingale and the limit theory of multipower variation devel- 
oped for semimartingales does not suffice to derive the desired kind of limit results. The 
basic tool we establish and apply for this is a normal central limit theorem for triangular 
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arrays of dependent Gaussian variables. As a case of some special interest for applica- 
tions, particularly in turbulence, the central limit behaviour of the realised variation 
ratio, that is, the ratio of bipower variation to quadratic variation, is briefly discussed. 
Some specific examples of choice of g are also considered. 

The turbulence context concerns time- wise observations of velocities at a single location 
x in space. More generally, it would be of interest to develop the theory of multipower 
variation corresponding to a setting where velocities are observed along a curve r in 
space-time. More specifically, suppose that velocity Y t (x) at position x and time t is 
defined by 

Y t (x)= f g(t-s t x-S)v B (£)W(dtd8), 

J A+(x,t) 

where W denotes white noise, <Jt(x) is a positive stationary random field on E 2 , j is a 
deterministic damping function and A is a subset of space-time involving only points with 
negative time coordinates. Then, with the curve r parametrized as t(w) = (x(w),t(w)), 
say, the problem is to study multipower variations of the process X defined as 

X w = f g(t(w)-s,x(w)-Z)a 8 (Z)W(dtds). 

J A+t(w) 

Among the questions that this raises is that of proper definition of filtrations. 

In another direction it would be of interest to extend the results of this paper to power 
and multipower variations of higher-order differences of Y. In particular, this might yield 
normal central limit theorems for the whole range of values of S and it could also lead to 
more robustness against drift processes. For some recent work on quadratic variation of 
higher-order differences, see [18, 19] and references given there. 

8. Proofs 

As our proofs are rather long and technical, let us briefly outline the scheme. The proofs 
basically consist of four main steps: 

(i) First, we provide a Wiener chaos decomposition for the functional V(p J 1 , . . . ,p 3 k )l l , 
j = 1, . . . ,d, that appears in Theorems 1 and 2 (see Section 8.3). 

(ii) In a second step, we prove the weak law of large numbers and the central limit 
theorem for the normalized version of (V(pj, . . . ,Pfc)™)i<j<d, using its Wiener chaos de- 
composition and recent techniques from Malliavin calculus derived in [30-32] (Sections 8.4 
and 8.5). 

(iii) In order to prove Theorems 3 and 4, we first show that the true increments A™y 
can be replaced by the quantity a^_iy n AfG without changing the asymptotic limits 
(see (8.18) and (8.26)). For this step, the conditions (4.6) and (4.8) on the measure 7r" 
are absolutely crucial. 
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(iv) In the last step, we apply the following blocking technique: We divide the interval 
[0,t] into big sub-blocks (whose lengths still converge to 0) and freeze the volatility 
process a at the beginning of each sub-block. Then Theorem 4 (respectively, Theorem 
3) follows from Theorem 2 (respectively, Theorem 1) applied to the family of functionals 
(V(G,p{, . . •,Pfe)")i<j<c2 and the properties of stable convergence. 

Below, all positive constants (which do not depend on n) are denoted by C, although 
they might change from line to line. 

8.1. Some elements of Malliavin calculus 

Before we proceed with the proofs of the main results, we review the basic concepts of the 
Wiener chaos expansion. Consider a complete probability space (f2, J 7 , P) and a subspace 
Hi of L 2 (il, J 7 , P) whose elements are zero-mean Gaussian random variables. Let H be 
a separable Hilbcrt space with scalar product denoted by (•, - )h and norm || • ||h- We will 
assume that there is an isometry 

h^W{h) 

in the sense that 

E[W(hi)W(h 2 )} = (hi,h 2 )m. 

It is easy to see that this map has to be linear. 

For any m > 2, we denote by H m the mth Wiener chaos, that is, the closed sub- 
space of i 2 (51,J-",P) generated by the random variables H m (X), where X £ Hi, 
E[X 2 } = 1 and H m is the mth Hermitc polynomial, that is, Hq[x) = 1 and H m {x) = 

We denote by 

I m : H 0m — > H m 

the isometry between the symmetric tensor product H 0m , equipped with the norm vm!|| ■ 
and the mth chaos H m ', see Section 1.1.1 in [29] for its definition. For h £ H® m , 
we set I m (h) := l m (h), where h is the symmetrization of h. For any g £ H® m , h £ H® n , 
n, m > 0, it holds that 

E[I m (g)I n (h)}=6 

HI II 

For any h = hi <8> ■ • • ® h m and g = gi (8 • ■ ■ ® g m & H 8 "™, we define the pth contraction of 
h and g, denoted by h® p g, as the element of H® 2 ( m_p ) given by 



h® p g= (hi,gi) a -- • (hp,g p }wh p+1 • • • ® h m ®g p+1 ® ■ ■ ■ ® g m . 
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This definition can be extended by linearity to any element of H® m . 

Now, let Q be the cr-field generated by the random variables {W(h)\h £ H}. Any 
square-integrable random variable F £ L 2 (Q,Q,P) has a unique chaos decomposition 

oo 

F=J2 I m(h m ), 

m=0 

where h m £ H 0m (see [29] for more details). 

Finally, we adapt the theory of Wiener chaos expansion to the set up of Section 3. 
Let Q be the cr-field generated by the random variables (^i,n)n>i,i<i<[nil an d T~L\ be 
the first Wiener chaos associated with (-?Q,n)n>i,i<i<[ntl; that is, the closed subspace of 
L 2 (il,Q,P) generated by the random variables pQ,n)n>i,i<i<[ni]- Notice that Hi can 
be seen as a separable Hilbert space with a scalar product induced by the covariance 
function of the process (-2Q,n)n>i,i<i<[ni]- This means we can apply the above theory of 
Wiener chaos expansion with the canonical Hilbert space M = Hi- Denote by H m the 
mth Wiener chaos associated with the triangular array (Xi : „)„> 11 < i <r 7l ti and by I m the 
corresponding linear isometry between the symmetric tensor product T-Lf m (equipped 
with the norm A/m7| ■ ||-%®">) an d the mth Wiener chaos. 

8.2. Preliminary results 

First of all, let us note that w.l.o.g. the volatility process a can be assumed to be bounded 
on compact intervals because a is cadlag and it is integrated with respect to W . This 
follows by a standard localization procedure presented in [6]. Furthermore, the process 
Ft, defined by (4.4), is continuous, because a is cadlag. Hence, F t is locally bounded and 
can be assumed to be bounded on compact intervals w.l.o.g. by the same localization 
procedure. 

Next we establish three lemmas. 

Lemma 1. Under assumption (4-4), it holds that 

E[\A^Y\p]<C p tP, i = 0,...,[nt], (8.1) 

for all p> 0. 

Proof. Recall that g' is non-increasing on (b, oo) for some b > 0. Assume w.l.o.g. that 
b > 1. Observe the decomposition 

KY= g(--s)a s W(ds)+ ^'Y g(- - s) - g ( — - s) )a s W(ds). 

J{i-l)/n \ n J J-oo V \ n J \ n J J 

Since a is bounded on compact intervals, we deduce by Burkholdcr's inequality 
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We immediately obtain the estimates 



ds < — . 



because g' is continuous on (0,oo) and a is bounded on compact intervals. On the other 
hand, since g' is non- increasing on (6,00), we get 

1 , ^ „f„\\ J2 F(i-l)/n 



+S ) ~- 9 ^) 4-1)/— «ds< 



The boundedness of the process F implies (8.1). □ 

Next, for any stochastic process / and any s > 0, we define the (possibly infinite) 
measure (recall that g(x) := for x < 0) 



K n , A) _ E 1a(9( x - IM -gW) 2 /Ld^ AeB(ffi). (8.2) 
Jo 00 ^ -1 /™)-^)) 2 ^ 



Lemma 2. Under assumption (4-4), it holds that 

sup < s ((e,oo))<^"((e,oo)) (8.3) 
«e[o,t] 

/or any e > 0, where ir n is given by (4-5). 

Proof. Recall again that \g'\ is non-increasing on (6,00) for some b > 0, and assume 
w.l.o.g. that b> e. Since the processes a and F are bounded we deduce exactly as in the 
previous proof that 



9\x ~g(x) o- s _ x dx 



x--j-g(x)j a 2 s _ x Ax+ f (g(x--j-g(x)j a 2 s _ x dx 



x-^)- g(x)\ dx + n 



This completes the proof of Lemma 2. □ 



Finally, we present the following technical lemma. 
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Lemma 3. Under the assumption (CLT), there exists a number I > 1 and positive se- 
quences e$ — > 0, j = 1, . . . , I, such that < < • • • < el o/ad 

4 1 )=o(n- 1 /( 2 7(fAi))) 5 ^((4 i ),oo))=o(7i- 1 /(f A1 )), (8.4) 

(4+ 1 )) 27 7r"(( £ « ) oo))=o(n- 1 /(f^)) ) j 1 / !. (8.5) 

where p = mhii<i<k,i<j<d(j4)- 

Proof. Assume first that p > 1. Recall that 7 > 1/2 . Set s„ — n _K J ' , j = 1 , . . . , I, with 1 > 
Ki > ■ ■ ■ > Ki > 0. The condition 7r™((e£f , 00)) = 0(n A ( 1-Kj ')) for some A < — 1, presented 
in (4.8), implies that conditions (8.4) and (8.5) are satisfied if we find 1 > Ki > ■ ■ ■ > Ki > 
such that 

KK1 + -, 

(1 + A) - KjX - 2k 3 - + i7 < 0, 1 < j < / - 1. 
From the first and the last inequality, we deduce by induction that 



1 1 + 2 7 V / 2 7 s 1 



must hold. 

When 27 > —A, the term on the right-hand side converges to 00 as / — > 00. In that case 
it is easy to find constants 1 > Ki > ■ ■ ■ > k; > such that (8.4) and (8.5) are satisfied. 

When 27 < —A, the limit of ^Etof - ! 1 )' i s x+2 ( as ' ~^ 00 ) anc ^ tnc restriction 
on Ki becomes 

1 1 + A 

7T < «i < 



2 7 A + 2 7 

Notice that ^ < j^^j because 7 > 1/2. The existence of the positive powers Kj, j = 
2, . . . ,1 that satisfy the original inequality follows by an induction argument. 
Assume now that p < 1 . Recall that 7 must satisfy 

1 

7> 7T 

2p 

and A < — -. Again the conditions (8.4) and (8.5) are satisfied if we find 1 > K\ > • • • > 
Ki > such that 

1 1 + Xp 

2~fp \p 

1 



A ) - KjA - 2k j+ i7 < 0, 1 < j < I — 1 
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Notice that the second inequality has solutions because A < — K Moreover, we deduce as 
above that the inequality 

^ <Ki< ^SW + {--) Ki 

must hold. Again the more complicated case is 2 7 < —A. By letting I — > oo, the restriction 
on Ki becomes 

1 1/p+A 
2 1P <Kl< A + 2 7 ' 

Note that 2^ < because 7 > ^ ■ As before, the existence of the positive powers 

Kj, j = 2, . . . , I that satisfy the original inequality follows by an induction argument. □ 



8.3. Some notation 

In this subsection we introduce various notation connected to the Wiener chaos expansion 
for the functionals V(p{, . . . ,p 3 k )™, j = 1, • • • , d, and present some first convergence results. 

Recall that the covariance matrix of (A !in , . . . , Xi+i >n ) is invertible for any I > 1 and 
n > 1. Let X™(1), . . . , A"(fc) be an i.i.d. ^(0,1) sequence that spans the same linear 
space as Xi <n , . . . ^X^k-i.n (such a sequence can be constructed by the Gram-Schmidt 
method). Thus, it has the representation 

k 

Xnj)=$>£X i+i _ 1>n) j = l,..., h (8.6) 
i=i 

for some real numbers a[* . Note that 

K\<c 

for all l,j,n, because E[Xf n ] = 1 for all i,n. 

For any l<j<d, we obtain the Wiener chaos representation 

oo / [nt] \ 

v( P {,. . . ,p{)? - P w ,t = J - - E + ( 8 - 7 ) 

Pn—>Pi z — ' V n * — ' J / 

m=2 \ i=l / 

where the fmj(i) G H 0m are given by 

fcie{i,...,fc} 
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for some coefficients fc (j). We set 

c(i)=n/^)iiH®".= E K,..,k m m 2 - (8-9) 

feie{i,...,fe} 

Note that 

OO 

var(|Xi,„^ • • • \X i+k _ 1>n \<) = J2 m!C(j) < C (8-10) 

m=2 

for all n, j, because E[Xf n ] = 1 for all i,n. Finally, when f^j(i), fcm (j) and c" n (j) 

correspond to some particular choice of powers pi,...,p k , we use the notation f^(i), 

c k u ...,k m and C- 

Now assume that the assumptions (3.8) and (3.9) of Theorem 2 hold. Since afj in 
(8.6) is a continuous function of r(l), . . . ,r n (k — 1) and the Gaussian process Q is non- 
degenerate, we have that 

and the sequence Qi(l), . . .,Qi(k) given by 

ftO) = y^%'Qw-l; J = 1, • - • , fe, (8.11) 
i=i 

is an i.i.d. iV(0, 1) sequence. Now, let us associate f m ,j(i), c fci,...,fe TO (i) an d c m (j) with 
the functional Vq(pJ, . . . - Prf pi t, where 

^....y—^IQiK 1 ---!^! 4 ]. 

by (8.7), (8.8) and (8.9). By a repeated application of the multiplication formula (see 
[29]), we know that cJJ fc (j) is a continuous function of r„ ( 1 ),..., r„(fc — 1). Since 
r n (j) -> P(i) we obtain 

^ 1 ,...,* m (j')->c fcl ,... lfcra (j), C(j)^c m (j), (8.12) 
(^B.i^(i + I)>H*« ^ + (8.13) 

covdX^J^ 1 • • • \X i+k _ ltn \^ , |X ijn |^ 2 • • • \X i+k _ ltn \pi 2 ) 

OO 

= E m H/^ 1 (l):/™ >J - 2 (l))H«"' 
m=2 

-> covdQ,^ 1 • • • \Q i+k -i\^ , \X hn f* ■ ■ ■ \Q i+k -i\ p *) 

OO 

= E m! (/m,ii(l),/m,j 2 (l))H®"'- 
m=2 



(8.14) 



Asymptotic results 

8.4. Proof of Theorems 1 and 3 



1183 



Proof of Theorems 1. Since V(pi, . . -,Pk)t is increasing in t and the process Pjn,...,p k t 

is continuous in t, it is sufficient to prove V(p±, . . . ,Pk)t ~ Ppi]—,Pk^ f° r a nxe( l t > 0. 
Note that 

I(/™(i),.C(i + 0)h^I<C 8 

icc(i), mi + o>h»»i < c^ m (Mor + • ■ ■ + ir„o - * + i)n, 

where the bounds are not comparable. Now, due to assumption (3.3), r(j) — ¥ as j — > oo . 
Thus, there exists an H such that ICr 1 / 2 ^' — k + 1)| < 1 for j > H (for any fixed C). By 
(8.15) we have (for any m > 2) 



E I + 0>H*» I < C + E K/m(l), + 0>H- I 

(=1 V l=H / 

(n-1 \ n-1 

H+Y,(C\r n (l)\) 2 <C?CE r W- 



(8.16) 



Hence, 



^ oo / n—1 \ 

vaxiVfa, . . . ,p fc )?) < - E ™ ! C 1 + E r ^ • 

m=2 \ (=1 / 

The latter converges to due to (8.10) and assumption (3.3). □ 

Proof of Theorem 3. Assume first that Z\ = 0. Recall that 

£[|A«yH<CV«, £[|A?Gf] < Crl (8.17) 

for any q > 0, due to Lemma 1. 

We assume for simplicity that k = \.p\ = p. The general case can be proved in a similar 
manner by (8.17) and an application of the Holder inequality. 

Since V(Y,p)" is increasing in t and the limit process is continuous in t. it suffices to 

prove the pointwise convergence V(Y,p)" — ^> (i p Jq |cr s | p ds. For any I <n, we have 
,t 1 [«*] 

nw-/%> / k 8 pd S =^E(i A " y i p -i a (-i)/» A ? G n +jR "' i > 

JO 



where 



. M [it] \ 

R t' l = ^ Ei^-i)/» a " g I p -EK-daI p E i a ? g i p 

=1 i-1 iC7,(i) / 



1184 



O.E. Barndorff-Nielsen, J.M. Corcuera and M. Podolskij 



and 



[It] 



[it] 



3=1 



+ M^r 1 Eko--i)/ i r-^ t Kl p d^ 



n V I 'I 



i>i- 



The assumption (LLN) implies that V(G,p)™ /x p t. Since er is cadlag and bounded on 
compact intervals, we deduce that 

lim lim P(\Rt' l \ >e) = 

for any e > 0. Hence, we are left to prove that 

["'] 



Ed A " y l p - k(i-u/«A?Gn -A o. 



By applying the inequality ||x| p — |y| p | < p\x — y\(\x\ p ^ 1 + for p > 1 and ||a;| p — 

\y\ p \ < 1^ — y\ p for P < 1: (8-17) and the Cauchy-Schwarz inequality, we can conclude that 
the above convergence follows from 



it] 



-L]>>[|A»y -<7 (i _ 1)/n A«G| 2 M0. 



1=1 



Observe the decomposition 



where 



5! 



g[ s (cr s - a {i _ 1)/n )W{ds), 

(i-l)/n \ n / 



r(*-i)/n 

0-(i-l)/ n / ( #( f -S ) -g( -3 ] ]W(ds), 



(i— l)/n— er 



c n,e = 



(i-l)/n- 



(8.18) 
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°~(i-l)/n 



(i—l)/n—e 



5 s 

n 



i-1 



s )W(ds). 



By Lemma 2 and the boundedness of a on compact intervals, we deduce 

[nt] 



and by (4.6) we obtain that 



[nt] 



lim J-Vi?[|Cn 2 ]=0. 



Next, we get 



'"n i=1 ""'n 



["*] r i/n 



E 



=1 -/(i-l)/n 



■ 9 s)(a s - o- (i _i)/„) ds 



Set f (s, 77) = sup-flay — ay | 2 |s, r G [— t, i], |r — s| < 77}. Then we obtain 



-1 [«*] -1 ["*] r / • 1 



Moreover, for any k > 0, since <r is cadlag, there exists n big enough that 
'i-1 



,n 1 I < K + (A(T(j_i)/ n ) l{(A<r (j _ 1)/n )2>K}, 



so 



[nt] 



[nt] 



^ E ^"l 2 ] < « + ^ E S[( A ( T( 4 -l)/„) 2 l { (A CT(l _ 1)/ „)^> K }] 

i=l 

- E ( Acr s) 2l {(A CT3 ) 2 >K} 



i=l 



<K + E 



then 



[nt] 



lim A-V£[K| 2 ]< 

i-livi nr2 u 1 1 J 



n->oo nr. 



and the convergence to zero follows, letting k tend to zero. 
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Finally, observe the decomposition B™' e = B™' e (l) + £?™' e (2) with 



bt (i) = i 9 ft " s ) - 9 (^r - s ) ) - *<">/~>*w. 

AT (2) = H- l)/n - e - ^ 1)/n ) f^_ s (g (1 - S ) " * -)) WW 



We obtain the inequalities 



it] 



J^£p^i)| 2 ]<i££ 

n i=l i=l 
[nt] [nt] 



i - 1 



i - 1 



1/2 



(8.22) 



By using the same arguments as above, we have that both terms converge to zero and 
we obtain (8.18), which completes the proof with Z\ = 0. 

To prove the general case, with Z\ ^ 0, we consider, for simplicity, the case k = 2. 

Assume first that < p t ,p 2 < 1. We have ||x x + yi\ Pl \x 2 + y 2 \ P2 - \yi\ Pl \y2\ P2 \ < 
C(\ Xl \ Pl \x 2 \ P2 + \xi\ Pl \y 2 \ P2 + \yi\ Pl \x 2 \ P2 ). Hence we deduce 



\V(Z, Pl ,p 2 )?-V(Y, Pl ,p 2 )?\ 



[nt]-l 



and the result follows by (4.10). 

Next, assume that p\ <p 2 , p 2 > 1. We deduce that 

|(T/(^Pi,P2)?) 1/p2 - (V(F,pi,p 2 )r) 1/p2 l 

^ C — E iA?^riAr + i^r + 



[nt]-l \ !/P2 

— ^ |A?Zir|A? +1 Z a |» j 



\ nr ™ i=i / 



1/P2 



which completes the proof of Theorem 3. 



□ 



Asymptotic results 

8.5. Proof of Theorems 2 and 5 
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Proof of Theorem 2. We first show the weak convergence of finite-dimensional distri- 
butions and then prove the tightness of the sequence ^Jn{V{p\, . . . , — , t)i<j<d- 
Step 1: Define the vector Z n (j) = (Z*(j) t Z c n {j)) T , 1 < j < d, by 



[nb,\ 

Z l n (j) = ^= {\Xi, n \^--\Xi +k -i,n\ vi -P^ J), (8-23) 

i=[nc ( ]+l 



where (q, &;], I = 1, . . . ,e, are disjoint intervals contained in [0,T]. Set Z l n = (Z l n (l) 7 . . . , 
Z l n (d)), I — 1, . . . , c. Clearly, it suffices to prove that 

(Z l n )i<i< c (P 1/2 (B bl - B Cl ))i<j< e > 
where the matrix (3 is given in Theorem 2. By (8.7), we have the representation 

[nb,] 



oo / [n&jj \ 

m=2 \ V i=fnc,l+l / 



Set F™ = -j= J2iL[nc,]+i /m,j(*)- By Theorem 2 in [5] we obtain the weak convergence 
of finite-dimensional distributions when we show that: 

(i) For any 1 < I < e, 1 < j < d, we have 

oo 

Jim limsup V m!||.F£ j(j')||h®* = 0. 

m— A/+1 

(ii) For any m > 2, 1 < Z < e and 1 < ji, ja < rf, we have constants Cjy such that 

lim m!(^,(ji),^,(j 2 )>H(8m = C m .i{j ll j 2 ) 1 

n— too 

and Em=2 Cm,l(h,j2) = Pj u h{ h l - c l) ■ 

(iii) For any 1 < ii 7^ Z2 < e and 1 < ji, j'2 < we have 

lim {FZ M {h),FZ M {h)) m rr, =0. 

(iv) For any m > 2, 1 < I < e, 1 < j < d and p = 1, . . . , m — 1 

lim K 1 iCj)®p^i(j')II|o (3 (»-,)=0. 
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Note that it is sufficient to prove (i), (ii) and (iv) for 1 = 1, h = l and a; = 0. In this 
case we use the notation F™(j) = 
(i) and (ii) As in (8.16), we have 

m!<K0i),-C(ja)>H*» 

= m\ ((f^ h (1), f£, h (1))h«- + \ ][> - 0(/^ (l),/:, 3! (1 + Ok®-) (8-24) 



n-l 



< 



Cm\{fZ, h (l),f^ h (l))m^ 1 + 



\ 2=1 / 

Since X)z*li r (0 < 00 ) we obtain by (8.12)-(8.14) and the dominated convergence theorem 
lim ml{F^tUi),K,lU2))B^ 

n— foo ' 

= C m (j 1 ,j 2 ) 

= ^(^,(1)^(1^ 

and X)m=2 Cm(ji;J2) = /8ji,j2 (notice that /8j lli3 \, is finite due to the dominated conver- 
gence theorem). Hence, we deduce (ii). On the other hand, we have 

oo oo 

limSUp E m! (/mj 1 ( 1 )Jm j2 ( 1 ))l» m = E m! (/mjl( 1 )Jmj 2 (l)) H ®'»<00. 
n ^°° m =N+l m=N+l 

Thus, we obtain (i) by (8.24). 

(iii) W.l.o.g. consider the case j = j\ = j 2 . For any h < l 2 , as in (8.16), we have 

[ra&jj [nb l2 ] 

\{K,iM^,iM^\<- E E 

/i=[nc( 1 ] + l i=[nc( 2 ] + l 

Assume w.l.o.g. that c\ 1 = 0, =ci 2 = l and fof 2 = 2 (the case ^ < ci 2 is much easier). 
Then, by condition (3.9), we obtain the approximation (as in (8.24)) 



/ n n— 1 \ 



3=1 3=1 



since E^U r 0') <°°- 

(iv) A straightforward computation shows that 

\\KU)®pKU)\\mHm-rt 
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n-l 

E - ialKfli* ~ *3|)C- p (Ki ~ i4|)C- p (|i 2 - js|) 



„ n-l 



n,i2,j3,i4 = l 



for some constant C mj . The latter is smaller than 



C ' 



E l<«IK(OIIC- p (K-ft|)IIC- p (l«-^l)l 



71 

i,h,l=l 



n—X / n— 1 



= -E Ei<«nc- p (K-^i)i 

h=l\i=l / 

Now, for any < e < 1, wc obtain by the Holder inequality 

n- 1 E ( E i<«iKr p (i^i)i) 2 

0</i<n-l ^0<t<n-l ' 

^ n_1 E ( E K(i)\K- p (\i-h\)\ 

Q<h<[ne] M)<i<n-1 

n-l /[ns/2] \2 

+ 2n_1 E E Kmrn- p {v-h\)\\ 

h=[ne] \ i=0 / 
n-l / n-l \ 

+ 2n ~ 1 E E i<( i )iic- p (N-^i)i 

fe=[ne] \h=[nc/2] J 



< c (4 E icwi) + E icwi) E icwi) 

^ \)<i<n-l ' \<i<n-l ' \ne/2]<h<n-l ' ' 

The latter is smaller (again by (3.9)) than 

/ / \ 2 / \ 2 P/ m / \ 2(m-p)/m\ 

4( E r «) +( E m ( E '(*)) ) 



that converges to C£(^^. r(«)) 2 as n — > oo. Thus, we obtain (iv) by letting e — > 0. 
Siep 2: Clearly, it suffices to consider the case d=l, p\ =pi- Set 



oo / [nt] \ 

V^(y(pi, . . . , Pk )? - p^...^) = E M ^ E + °p(^ 1/2 ) 

m=2 \ * n i=l / 

= :^ l + 0(n- 1 / 2 ) 
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(where the approximation holds locally uniformly in t) and 



N / 1 [nt] > 

^ = EM-iE/™(*) 

m=2 \ V i=l y 



In step 1, we have proved that conditions (i)-(iii) of Theorem 2 in [5] arc satisfied. 
Then by (2.3) of Theorem 2 in [5] and the Cauchy-Schwarz inequality, we obtain the 
approximation 



< 



A 4 



/3 2 ([nt] - [r^Mnfa] - [nt]) /3 2 (t 2 - ^) 2 
" ^ " A 4 

for any t\ <t<tz and A > 0. On the other hand, (8.14) and (8.24) imply that 

lim E[\Z?-Z?< N \*] = 

N—foo 

for any n and any t. Using this we conclude that 



p{\z™-zi\>\\zt-z*\>\)<c f - 



A 4 



for any ii < i < t 2 and A > 0, from which we deduce the tightness of the sequence by 
Theorem 15.6 in [20]. This completes the proof of Theorem 2. □ 

Proof of Theorem 4. We only consider the case Z\ = (if Z\ ^ 0, we can proceed as 
in Theorem 3). Also, for the sake of simplicity, we take d=l, k = 1 and pi =p. We use 
the decomposition from the proof of Theorem 3: 

I [It] [It] \ 

= M^/2\°U-D/i\ P E MG\r-» p r^Wu-W\ P ) ( 8 - 25 ) 
V ™ 3=1 iehU) i= x ' 

[nt] 

2 — 1 



for any / < n, with 

[it] 



$ ,t = -k* Ei^-D/« A " G r-Ei%-D/'i p E i a " g i p 

VnTn \ 4 =i 3=i iGJ,(i) y 
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[It] 



/ [«] ,t 
\/«Mph" i y'|a'(j-i)/i| P - / k« 

V ,-=i Jo 



Observe that under the assumption (CLT) we obtain the weak convergence 

V^(V(G, P )?~fi p t)A^pB t 

(see Theorem 2). Since E[Gt(V(G,p)™ — n P t)\ = for any t > 0, because G has a sym- 
metric distribution, we deduce (by Theorem 5 in [3]) that 

(G t ,VK(V(G,p)? - M p*)) 4 (G u yfpB t ). 
Now, an application of the condition D" from Proposition 2 in [1] shows that 

VE(V(G,p^-fi p t)^^B t . 
By the properties of stable convergence, it follows immediately that 
/ 1 M [it] \ _ [it] 

V " j=i ieijO') i =1 ' J =1 

for any fixed I. On the other hand, we have 

[it] 



.7=1 Jo 



as Z — > oo . 

Now we need to prove that the other summands in the decomposition (8.25) are neg- 

n i 

ligible. The negligibility of the term R t ' is shown as in the proof of Theorem 7 in [3] 
but by using condition (4.7) instead of the Holder continuity of index 7. So we are left 
to prove that 

[nt] 

E(I A ?>T - k (i _ 1)/n A?Gf ) -A 0. (8.26) 

v 1—1 

By applying, for p > 1, the inequality \\x\ p — \y\ p \ < p\x — y|(|a;| p_1 + |y| p_1 ), (8.17) and 
the Cauchy-Schwarz inequality, and, for p < 1, — < \x — and the Jensen 
inequality, we have 

[nt] [nt] 

^siiAryr-k (i _ 1)/n A«Gn<— w £(£|A?F-a (i _ 1)/n A?q 2 )<^/ 2 . 

v 1=1 V rt j = l 
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Now we use a similar decomposition as presented in the proof of Theorem 3: 

I 

A?Y - £7(i_i)/ n A"G = ^ + Si ' " + E ^ " ' " , 

3=1 

where A", _B™' e ™ are defined as above, < < ■■■ < £n < ei' +1 ' = oo and 

(i-l)/ra-e<?> 



l(i-l)/n-s 

By the assumption (CLT) and Lemma 2, we obtain the following inequalities (since a is 
bounded on compact intervals) 



[nt] 



l 7 (pAl) + l/2 



iT n . =1 
[nt] 

i^E(^i B r £ " , i 2 ) (pA1)/2 <^ 1/2 i4 l) r (pA1 



[nt] 

-t^mE^I^"^^ ?+1 Y)^ 1)/2 <cnV^|4^) r ^ 1 )| vr n ([ 4,o i 4^) ) |^i)/ 2) 

j = l,. ..,1-1, 

[nt] 

— L^E^I^ 4 ' 4 l 2 )^^ 2 < Cn^((40 jOo)) (^)/ 2 . 
V nr « i= i 

Then we deduce by (CLT) and Lemma 3 

[nt] 



-i^El^-^-i)/^"^^ - 

v i— 1 



This completes the proof of Theorem 4. □ 
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